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ABSTRACT 

The massive OB-type binary a Ori AB is in the centre of the very young a Orionis 
cluster. I have computed the most probable distances and masses of the binary for sev- 
eral ages using a dynamical parallax-like method. It incorporates the BVRIH-ha,nd 
apparent magnitudes of both components, precise orbital parameters, interstellar ex- 
tinction and a widely used grid of stellar models from the literature, the Kepler's third 
law and a minimisation. The derived distance is 334I22 PC for an age of 3±2 Ma; 
larger ages and distances are unlikely. The masses of the primary and the secondary 
lie on the approximate intervals 16-20 and 10-12 M©, respectively. I also discuss the 
possibility of a Ori AB being a triple system at ~385pc. These results will help to 
constrain the properties of young stars and substellar objects in the a Orionis cluster. 
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1 INTRODUCTION 

The Trapezium-like system a Ori, that illuminates the en- 
colure of the Horsehead Nebula, is the fourth brightest star 
in the young Ori OB 1 b association. The multiple system 
is composed of at least five early-type stars (Burnham 1892; 
Greenstein & Wallerstein 1958; van Loon & Oliveira 2003; 
Caballero 2007b). The two hottest components, a Ori A and 
B (09.5V and B0.5V), are separated by only ~0.25arcsec 
and were for a long time "the most massive visual binary 
known" {Ma+Mb ~ 25 + 15 Mq- Heintz 1974). Although 
the binary has not yet completed a whole revolution, the 
orbital parameters are relatively well determined (Hartkopf, 
Mason & McCalister 1996; Heintz 1997; Horch et al. 2002). 
It has been suggested that a Ori AB is a hyerarchical triple, 
being the primary a short-period, double-line spectroscopic 
binary (Frost & Adams 1904; Henroteau 1921; Miczaika 
1950; Bolton 1974; Morrell & Levato 1991). However, a large 
amount of accurate, comprehensive spectroscopic investiga- 
tions have failed to confirm this hypothesis (Heard 1949; 
Conti & Leep 1974; Humphreys 1978; Bohannan & Gar- 
many 1978; Garmany, Conti & Massey 1980; Simon-Diaz & 
Lennon, priv. comm.). 

The a Ori system is located in the centre of the well- 
known a Orionis open cluster. The proper motions, radial 
velocities and spacial distribution of stars in this cluster 
strongly suggests a physical association between a Ori it- 
self and the young cluster (Zapatero Osorio et al. 2002a; 
Caballero 2007a, 2007c - see also Jefi'ries et al. [2006], who 
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discovered a second older and kinematically and spacially 
distinct population). Because of its youth, comparative near- 
ness and low extinction, the cluster has become the rich- 
est hunting ground for brown dwarfs and planetary-mass 
objects in the whole sky (Bejar et al. 1999; Zapatero Os- 
orio et al. 2000, 2002b, 2007). There is plentiful material 
in the literature about this cluster, covering topics like the 
initial mass function down to a few Jupiter masses (Bejar 
et al. 2001; Gonzalez-Garci'a et al. 2006; Caballero et al. 
2007), jets and Herbig-Haro objects (Reipurth et al. 1998; 
Andrews et al. 2004), the frequency of accretors and discs 
(Zapatero Osorio et al. 2002a; Oliveira, Jeffries & van Loon 
2004; Kenyon et al. 2005; Oliveira et al. 2006; Hernandez 
et al. 2007; Caballero et al. 2007), the X-ray emission from 
young objects (Walter et al. 1997; Sanz-Forcada et al. 2004; 
Franciosini, Pallavicini & Sanz-Forcada 2006) or their pho- 
tometric variability (Caballero et al. 2004; Scholz & Eisloffel 
2004) . The most used values of heliocentric distance and age 
of the a Orionis cluster are d ~ 360 pc and ~3 Ma (Brown, 
de Geus & de Zeeuw 1994; Ferryman et al. 1997; Zapatero 
Osorio et al. 2002a; Oliveira et al. 2002) . There is a consen- 
sus in the literature that the cluster is younger than 8 Ma 
and older than 1 Ma. There is, however, a strong divergence 
of opinion on the heliocentric distance. Caballero (2007a) 
compiled determinations in the literature of the distance to 
the a Orionis cluster from the 352]Iijgg pc from Hipparcos 
parallax to almost 500 pc from colour-magnitude diagrams. 

Apart from the uncertainties of theoretical isochrones 
at very young ages, the derivation of the initial mass func- 
tion of the cluster is strongly affected by the uncertainty 
in the actual age and heliocentric distance (Jeffries et al. 
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Table 1. Photometry of cr Ori A and B from the literature. 



Band 




CT Ori A 


(7 Ori B 


References 


B 


[mag] 


3.85±0.05 


5.18±0.05 


Ca07a 


V 


[mag] 


4.10±0.03 


5.34±0.10 


tBrOO 


R 


[mag] 


4.15±0.04 


5.49±0.13 


tBrOO 


I 


[mag] 


4.41±0.04 


5.66±0.16 


tBrOO 


H 


[mag] 


4.81±0.10 


6.02±0.10 


Ca06 



2006; Caballero et al. 2007). There are other investigations 
that require a precise age determination, such as the evo- 
lution of the angular momentum due to discs and stellar 
winds (EislofFel & Scholz 2007), disc dissipation (Hernandez 
et al. 2007) and evolution of hot massive stars. In this Let- 
ter, I revisit a well known method for distance determina- 
tion: the dynamical parallax (e.g. Russell 1928). I apply it 
to the a Ori AB binary using state-of-the-art data and tools 
to determine its mass, age and heliocentric distance. 



2 ANALYSIS AND RESULTS 

The determination of the dynamical parallax of a binary 
of known orbital period, P, and angular semimajor axis, a, 
uses the Kepler's third law and a power-law mass-luminosity 
relation (e.g. Reed 1984). In this Letter, I instead use the 
grids of theoretical models from the Geneva group (Schaller 
et al. 1992). On the contray to other widely used grids, 
like those by the Lyon and Padova groups (Baraffe et al. 
1998; Girardi et al. 2000), the Geneva grids tabulate abso- 
lute magnitudes in a large number of passbands and ages 
down to 1 Ma and are valid up to very high masses (i.e. 
Ai > 10 Mq). Although there are more binaries and bi- 
nary candidates in the cluster (Caballero 2005; Kenyon 
et al. 2005; Caballero et al. 2006 and references therein), 
cr Ori AB is the only pair whose orbital parameters are 
known. Here, I have employed the parameters given by 
Hartkopf et al. (1996). The almost face-on orbit is charac- 
terised by a long period (P — 155.3±7.5a), a close angular 
separation (a = 0.2642±0.0052 arcsec) and a low eccentric- 
ity (e = 0.051±0.015). The orbital parameters are consistent 
with those of Heintz (1997; P = 158 a, a — 0.265 arcsec, e 
= 0.06). 

Using the standard units Mq, a (annum) and AU for 
Ma + Mb, P and the physical semimajor axis a, respec- 
tively, the Kepler's third law takes the simple expression 
(A^A + A4b)P^ ~ a^. Accounting for a — dtana ~ da, 
where d is the heliocentric distance in parsecs, and replac- 
ing the values of P and a from Hartkopf et al. (1996), then 
the third law for a Ori AB can be written as: 



Ma+Mb = 7.45 10" 'd'' 



(1) 



(A^totai = A^a+A^b). At the Hipparcos distance d = 352 pc, 
the total mass of the binary would be about 32 Mq , that 
is less than the classical value of total ^ 40 Mq, but is 
consistent with the value of Al total ~ SOA^q estimated by 
Caballero (2007a). 

I tabulate in Table [1] the BVRIH-hand magnitudes 
of both components in a Ori AB, taken from the litera- 
ture (ten Brummelaar et al. 2000 -tBrOO-; Caballero 2006 
-Ca06-; Caballero 2007a -Ca07a-). Except for the B-band 



measurement, that was estimated from the Tycho-2 BtVt 
magnitudes, all the data come from adaptive optics obser- 
vations. The Hipparcos catalogue also tabulated the non- 
standard Hp-hand magnitudes. The magnitudes and colours 
of both stars correspond to what it was expected of an 09.5V 
and a B0.5V at d ~ 350 pc. 

I have computed through a simple minimisation method 
which are the most probable heliocentric distances for sev- 
eral cluster ages. In particular, I have looked for the minima 
of the following chi-square distributions: 



x\d,MA,MB)^xl{d,MA)+xUd,MB) 



(for fixed ages and metallicities) , where Xa Xb s^re: 
^ ("^A,B - m^^e)^ 



2 ("iA,A - a) 2 \ - 

>^A=Z. 5^. ' ^B = 2^ 



'A, A 



5m\ 



(2) 



(3) 



(A = B, V, R, I, H). mx^[A,B] and 5mx,iA,B] are the 
observed apparent magnitudes and corresponding uncer- 
tainties of a Ori A and B in Table [TJ and gj = 
[A B] (i^i -'^[A.B]) are the theoretical apparent magnitudes 
that a hypothetical star of mass 7Vf[A,B] would have at an 
heliocentric distance d. To compute [A b] ' ^ have used: 
(i) the theoretical absolute magnitudes A'Ix,[a,b] from the 
basic grids of non-rotating stellar models with solar metal- 
licity (Z = 0.020), overshooting and OPAL opacities of the 
Geneva group, (ii) the colour excess E{B — V) = 0.05 mag of 
a Ori AB from Lee (1968) and (Hi) the interstellar extinc- 
tion law parameters A\/Av and Rv from Rieke & Lebofski 
(1985). In detail, I have used the grids with standard mass 
loss M and ages 1.0, 2.0, 3.2, 4.9 and 10.0 Ma (Schaller et al. 
1992 - Sc92) and high mass loss 2 x Af and age 3.2 Ma 
(Meynet et al. 1994 - Me94). The models do not provide 
data for other ages less than 32 Ma. Caballero (2006) mea- 
sured an average solar metallicity of solar-like stars in the 
cr Orionis cluster, [Fe/H] = O.OiO.ldex, which justifies the 
use oi Z = 0.020. 

For a fixed age (and metallicity), there is a one-to-one 
correspondence in the grids between M;^ [a,b] and A1[a,b] 
if both A and B stars are in the main sequence. For each 
heliocentric distance, there is only one corresponding to- 
tal mass of the hypothetical binary, A^ total <x dP (Eq. [1]). 
If the mass of the primary, A^a, is fixed, then the mass 
of the secondary is obtained from the simple expression 
A^B(d, A^a) ~ Al total (d) — A^A- To sum up, there is a value 
of for each trio (d, A^a, age). Because of the known spec- 
tral types of cr Ori A and B, I conservatively imposed that 
the masses of the primary and of the secondary should lie on 
the intervals 10 ^ Ma sS 50 Mq and 1.5 Af© < Mb ^ 
AIa, respectively. These constraints speed up the minimiza- 
tion but do not influence the results (see below). 

Fig. H] illustrates the x^ minimisation for an age of 
3.2 Ma and standard mass loss. Top and bottom windows 
display different coverage resolutions of the (d, A^a) plane. 
The results for normal and high mass loss for 3.2 Ma are 
almost identical. The plots for the other ages except for 
10.0 Ma are quite similar, with the minimum shifted along 
the X axis (distance). There are no solutions (i.e. the masses 
of A and B simultaneouly satisfy Eq. [T] and the above con- 
straints) for d < 250 pc and d > 500 pc. In addition, the 
4.9 Ma models do not provide solutions for d > 370 pc. 
Finally, there are no solutions at all for 10.0 Ma. At this 
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Figure 1. heliocentric distance for several masses of 

(7 Ori A. The dotted lines indicate curves of constant mass. The 
basic grid of models is from Schaller et al. (1992), the age is 3.2 Ma 
and the metallicity is Z = 0.020. Top window: whole interval of 
heliocentric distance in steps of 1 pc from 200 to 550 pc and of 
mass of the primary in steps of 2 Mq from 10 to 50 Mq . The 
minimum is for A^a = 18 Mq {d = 336 pc). Bottom window: 
same as top window, but for intervals of heliocentric distance in 
steps of 0.2 pc and of mass of the primary in steps of 0.2 Mq ■ The 
minimum in x^ is better constrained (cf. Table [2}. 



age, the primary has left the main sequence and got much 
brighter than the secondary. 

The values of {d, Ma) that minimise for given age 
and mass loss are provided in Table [2] The uncertainties in 
the values of d account for the error bars in the photometric 
data in Table[T]and in the orbital parameters P (5 %) and a 
(2%) in Hartkopf et al. (1996), and the size of the steps in 
the high resolution minimisation (Ad = 0.2 pc). The uncer- 
tainty in the masses is set to the step size, AA^a = 0.2 Mq. 
The results would be identical if no mass constraints were 
set. For an age interval of 3±2 Ma, the corresponding helio- 
centric distance interval is rf = 334I22 PC Distances larger 
than 400 pc and less than 290 pc are less likely for the 1-5 Ma 
age range; distances larger than 450 pc are highly unlikely 
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Table 2. Best fits of the dynamical parallax of a Ori AB. 



Mass 


Age 


d 


Ma 


Mb 




loss 


[Ma] 


[pc] 


[Mq] 


[Me] 




M 


1.0 


346±13 


20.1±0.2 


11.7±0.2 


8.70 


M 


2.0 


337±13 


18.3±0.2 


11.1±0.2 


9.38 


M 


3.2 


334±13 


17.6±0.2 


10.8±0.2 


9.76 


M 


4.9 


325±13 


16.0±0.2 


10.2±0.2 


10.56 


M 


10.0 










2x M 


3.2 


333±13 


17.4±0.2 


10.8±0.2 


9.72 



for all ages. Finally, I have determined the most probable 
masses of A and B for the Hipparcos parallax distance (d = 
352pc): Ma = 21AMq, Mb = 12.0 Mq for an age of 1 Ma. 
The minimum for 1 Ma is an order of magnitude smaller 
than for the other ages tested, meaning that a binary age 
older than 1 Ma would be unlikely at the Hipparcos distance. 



3 DISCUSSION 

The theoretical effective temperatures that correspond to 
the optimal fits lie on the intervals T^g = 30.4-34.6 kK for 
the primary and Tcff — 25.2-27.5 kK for the secondary. The 
hottest temperatures are for the youngest ages. These values 
are consistent with the expected Tcff for 09.5V and B0.5V 
stars, respectively (e.g. 09.5V: 30-35 kK - Popper 1980; Gu- 
lati, Malagnini & Morossi 1989; Castelli 1991; Vacca, Gar- 
many & ShuU 1996; Martins, Schaerer & Hillier 2005), and 
with previous measurements of the T^s of a Ori A (30.0- 
33.0 kK - Morrison 1975; Underhill et al. 1979; Morossi & 
Crivellari 1980; Repolust et al. 2005). The corresponding 
theoretical gravities (logg ~ 4.00-4.22) are also normal for 
class V at such temperatures. 

The minima of x^ in Table [5] are very sensitive to the 
variations of heliocentric distance and of mass: on the one 
hand, at fixed mass and age, a fluctuation of d of barely 
30 pc results in a change of three orders of magnitude in x^; 
on the other hand, at fixed distance and age, a fiuctuation of 
Ma of barely 5 Mq results in a change of almost two orders 
of magnitude in x^. The minima of x^ ^-re, however, quite 
unresponsive to the variations of the age between 1.0 and 
4.9 Ma (see last column in Table [2]). A younger age gives a 
slightly better fit results and that favours a slightly larger 
distance (d ~ 350 pc). The results do not strongly suggest 
a younger age of 1 Ma, but they are useful in excluding an 
older age. The absence of a solution at 10 Ma agrees with 
previous upper limits on the ages of the Ori OB 1 b as- 
sociation from the presence of very early-type stars in the 
main sequence (Blaauw 1964) and of the a Orionis cluster 
from spectral synthesis surrounding the Li I A670.78 nm line 
(Zapatero Osorio et al. 2002a). 

The derived distance interval for 3±2 Ma, d = 
334^22 P^i is consistent with the canonical distance to the 
a Orionis cluster of d ~ 360 pc, but is difficult to conciliate 
with the distance of 440 pc for 2.5 Ma that Sherry, Walter & 
Wolk (2004) used. The derived distance interval also devi- 
ates from very recent determinations of the distance to some 
elements in the Ori OB 1 complex. In particular, Terrell, 
Munari & Siviero (2007), using the eclipsing spectroscopic 
binary VV Ori close to Alnilam (e Ori) in Ori OB 1 b, and 
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Sandstrom et al. (2007), employing the Very Large Baseline 
Array in the Orion Nebula Cluster in Ori OB 1 a, have deter- 
mined very accurate heliocentric distances of 388-389 pc (see 
also Menten et al. 2007). These values are also lower than 
the classical distance to the Ori OB 1 complex of ~440 pc 
from average Hipparcos parallax (Brown, Walter & Blaauw 
1999; de Zeeuw et al. 1999). Because of projection effects 
and the large physical size of Ori OB 1 (Reynolds & Ogden 
1979), a Ori could be easily contained within the complex. 
Given their kinematic and spacial association, the a Ori 
system and the young a Orionis cluster are likely at the 
same heliocentric distance and also age, if one assumes that 
massive and low mass star formation in a cluster is coeval 
(Prosser et al. 1994; Massey, Johnson & Degioia-Eastwood 
1995; Stauffer et al. 1997; see, however, Sacco et al. 2007). 
Recently, it has been suggested that the a Orionis cluster 
is actually kinematically distinct from the Ori OB 1 b asso- 
ciation (Jeffries et al. 2006), just as 25 Ori is distinct from 
Ori OB 1 a (Bricefio et al. 2007). 

If a Ori AB were a hyerarchical triple, as described 
in Section [T] it would be located at a larger heliocen- 
tric distance. The hypothetical companion to a Ori A, to 
which I tentatively call a Ori F, would be 0.5 mag fainter 
than the primary in the 370-493 nm interval according to 
Bolton (1974). This wavelength interval corresponds to the 
U, B Johnson bands. Taking into account jtia = ttia+f -|- 
2.5 log ^1 -I- 10 *2.5 and the difference -Ba+f — Bb in 
Table [T] then the apparent magnitudes in the blue band 
of the three components would be related to the combined 
magnitude iJA+p through: _Ba ~ -Ba+f + 0.53 mag, Bb. ~ 
Ba+f + 1-33 mag and -Bf ~ Ba+f + 1-03 mag. I have looked 
for the distances and theoretical masses whose correspond- 
ing apparent magnitudes match the -B[a,b.fi relations and 
the Kepler's third law, {Ma + Mf) + Mb = 7.45 lO^'^d^ 
There are only a few solutions that simultaneously verify 
Tcff.A ~ 30.0-33.0 kK and Teff ,b « 26.0-30.0 kK, as expected 
from the spectral types of a Ori A-l-F and a Ori B. In the 
triple scenario, the F component would have an intermedi- 
ate temperature between the A and B stars (roughly BO.OV) 
and would orbit very close to a Ori A. The valid solutions 
are found for the narrow distance interval 370-400 pc and 
only for ages between 1.0 and 4.9 Ma. Although distances 
less than 290 pc and larger than 450 pc are ruled out again, 
the most probable distance to a Ori under the triple hy- 
pothesis, d ~ 385 pc, agrees very well with those of VV Ori 
and the Orion Nebula Cluster. Further high-resolution spec- 
troscopic studies are needed to ascertain the existence and 
characteristics of a Ori F. 



4 SUMMARY 

I have determined the most probable heliocentric distances 
and masses of the components in the young binary cr Ori AB. 
The used methodology is an improvement of the dynamical 
parallax method using a minimisation of observed and 
theoretical apparent magnitudes of A and B. The values 
range in the intervals 325 pc < d < 346 pc, 16 Mq < Ma ^ 
20 Mo, 10 Mq < Mb < 12 M© for 1-5 Ma. Ages and dis- 
tances larger than or equal to 10 Ma and ^^450 pc are ex- 
cluded from the minimisation. The theoretical effective tem- 



peratures of both components are consistent with the ob- 
served spectral types. Accounting for uncertainties in the 
orbital parameters and photometric data from the litera- 
ture, the derived distance interval for age — 3±2 Ma is d = 
334^22 pc. The values of are minimum for the largest dis- 
tances within this interval, that translate into the youngest 
ages. If there were a third component, a Ori F, in a tight 
orbit to a Ori A, then the system could be at a larger he- 
liocentric distance of about 385 pc. 

The a Ori star system is in the centre of the young 
a Orionis cluster. The knowledge of the age and heliocentric 
distance of the cluster is fundamental for the study of the 
initial mass function down to the planetary regime and the 
evolution of discs and angular momentum in very young ob- 
jects. 
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